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Abstract. In this paper, the notion of p-adic multiresolution analysis 
(MR A) is introduced. We use a "natural" refinement equation whose so- 
lution (a refinable function) is the characteristic function of the unit disc. 
This equation reflects the fact that the characteristic function of the unit 
disc is the sum of p characteristic functions of disjoint discs of radius p _1 . 
The case p — 2 is studied in detail. Our MRA is a 2-adic analog of the real 
Haar MRA. But in contrast to the real setting, the refinable function gen- 
crating our Haar MRA is periodic with period 1, which never holds for real 
refinable functions. This fact implies that there exist infinity many different 
2-adic orthonormal wavelet bases in £ 2 (Q2) generated by the same Haar 
MRA. All of these bases are constructed. Since p-adic pseudo-differential 
operators are closely related to wavelet-type bases, our bases can be inten- 
sively used for applications. 



1. Introduction 

1.1. p-Adic wavelets and pseudo-differential operators. According to 
the well-known Ostrovsky theorem, any nontrivial valuation on the field Q is 
equivalent either to the real valuation \ ■ \ or to one of the p-adic valuations 
| • | p . We recall that the field Q p of p-adic numbers is defined as the completion 
of the field of rational numbers Q with respect to the non- Archimedean p-adic 
norm | • \ p . This norm is defined as follows: if an arbitrary rational number 
x 7^ is represented as x = p 72 ^, where j = j(x) G Z, and m and n are not 
divisible by p, then 

(1.1) \x\p=P~\ x^O, \0\ P = 0. 

This norm in Q p satisfies the strong triangle inequality |x+?/|p < max(|a;|p, \y\ p ). 

Thus there are two equal in rights universes: the real universes and the 
p-adic one. The latter has a specific and unusual properties. Nevertheless, 
there are a lot of papers where different applications of p-adic analysis to 
physical problems, stochastics, cognitive sciences and psychology are stud- 
ied [H] [ID], P2] [IS], [21] [22] ( see a l so t ne references therein). In view 
of the Ostrovsky theorem such investigations not only have great interest in 
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itself, but lead to applications and better understanding of similar problems 
in usual mathematical physics. 

We recall that there exists a p-adic analysis connected with the mapping Q p 
into Q p and an analysis connected with the mapping Q p into the field of com- 
plex numbers C, there exist two types of p-adic physics models. For the p-adic 
analysis related to the mapping Q p — > C the operation of partial differentia- 
tion is not defined, and as a result, large number of models connected with 
p-adic differential equations use pseudo-differential operators and the theory 
of p-adic distributions (generalized functions) (see the above mentioned pa- 
pers and books). In particular, fractional operators D a are extensively used 
in applications (see fore-quoted papers and especially [M]). 

It is well known that the theory of p-adic pseudo-differential operators (in 
particular, fractional operators) and equations closely related to wavelet type 
bases. It is typical that p-adic compactly supported wavelets are eigenfunc- 
tions of p-adic pseudo-differential operators [3] [5], [16], [17], [18], [20] - [22]. 
Thus the wavelet theory plays a key role in application of p-adic analysis and 
gives a new powerful technique for solving p-adic problems. This theory starts 
development only in resent years and has many open problems. 

In [20], S. V. Kozyrev constructed the orthonormal compactly supported 
p-adic wavelet basis fll.2p in £ 2 (Q P ): 

(1.2) B lja {x) = p- l/2 X P (p- l j(p^x - a))n(\p^x - a\ p ), x G Q p , 

j G J p = {1, 2, . . . ,p — 1}, 7 G Z, a G I p = Q p /Z p . Kozyrev's wavelets (II. 2p are 
eigenf unctions of the Vladimirov fractional operator [3H IX]. Further develop- 
ment and generalization of the theory of such type wavelets can be found in 
the papers by S. V. Kozyrev [21], [22], A. Yu. Khrennikov, and S. V. Kozyrev 
[TB] , [PT] . J. J. Benedetto, and R. L. Benedetto [8], and R. L. Benedetto [9]. 
In [3], the multidimensional p-adic wavelets generated by direct product of 
the Kozyrev one-dimensional wavelets were introduced. In [18] , a new type of 
p-adic multidimensional wavelet basis was introduced: 

e^l{x) = p^ /2 x P {s{p^x - a))Q(\p^x - a\ p ), x G Q p , 

where s G J p;m , 7 G Z, a G I p . Here J p . m = {s = p~ m (s + Sip + • • • + 
Sm-ip" 1 ' 1 ) ■ sj = 0, l,...,p - 1; j = 0, l,...,m - l;s ^ 0}, m > 1 is 
a fixed positive integer. The multidimensional wavelets from [3] are a par- 
ticular case of the last wavelets. Moreover, in [3], [18], there were derived 
the necessary and sufficient conditions for a class of multidimensional p-adic 
pseudo-differential operators (including fractional operator) to have such mul- 
tidimensional wavelets as eigenfunctions. 

It remains to point out that for pseudo-differential operators from [3], [18] 
a "natural" definition domain is the Lizorkin spaces of distributions $'(Q™), 
introduced in [3J. The space $'(Q") is invariant under the mentioned above 
pseudo-differential operators. Moreover, the above mentioned p-adic wavelets 
belong to the Lizorkin space $(Q") of test functions. Recall that the usual 
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Lizorkin spaces were studied in the excellent papers of P. I. Lizorkin [24] , [25] 
(see also [23], |3Q|). 

It's interesting to compare appearing first wavelets in p-adic analysis with 
the history of the wavelet theory in real analysis. In 1910 Haar [12] constructed 
an orthogonal basis for £ 2 (IR) consisting of the dyadic shifts and scales of one 
piecewise constant function. A lot of mathematicians actively studied Haar 
basis, different kinds of generalizations were introduced, but during almost the 
whole century nobody could find another wavelet function (a function whose 
shifts and scales form an orthogonal basis). Only in early nineties a method 
for construction of wavelet functions appeared. This method is based on the 
notion of multiresolution analysis (MRA in the sequel) introduced by Y. Meyer 
and S. Mallat [28], [25], [27]. Smooth compactly supported wavelet functions 
were found in this way, which has been very important for some engineering 
applications. In this paper we introduce MRA in £ 2 (Q P ) and present a concrete 
MRA for p = 2 being an analog of Haar MRA in Z^(R). The same scheme as 
in the real setting leads to a Haar basis. It turned out that this Haar basis 
coincides with Kozyrev's wavelet system. However, 2-adic Haar MRA is not 
an identical copy of its real analog. In contrast to Haar MRA in £ 2 (M), we 
proved that there exist infinity many different Haar orthogonal bases in £ 2 (Q 2 ) 
generated by the same MRA. 

1.2. Contents of the paper. In Sec. |2J, we recall some facts from the p-adic 
theory of distributions [IT] . [32] . [33] . [34] . In Sec. El some facts from the 
theory of the p-adic Lizorkin spaces [3] are recalled. 

In Sec. HI by Definition 14.11 we introduce the MRA adapted to the p-adic 
case. In Subsec. I4.2[ we introduce the refinement equation (14.71) 

P_1 1 r 

whose solution <f){x) = f2(|x| p ) is the characteristic function of the unit disc, 
where where Q(t) is the characteristic function of the interval [0, 1]. The con- 
jecture to use the above equation as the refinement equation was proposed 
in [18]. The above refinement equation is natural and reflects the fact that the 
characteristic function of the unit disc B is represented as a sum of p 

pieces characteristic functions of the disjoint discs £>_i(r), r = 0, 1, . . . ,p — 1 
(see (EZD). 

In Subsec. I4.3[ the 2-adic MRA is constructed. Namely, we proved that 
MRA is generated by a refinable function which is the characteristic function 
4>{x) = Q[\x\ 2 ) of the unit disc B = {x : \x\ 2 < 1} C Q 2 and satisfies the 
refinement equation (14. 8 p 

1 1 r 
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By our MRA we construct 2-adic orthonormal wavelet basis (14.151) in £ 2 (Q 2 ), 
which is the Kozyrev basis (II. 2p for the case p = 2. It turned out that the 
Kozyrev wavelet basis is not unique orthonormal wavelet basis. 

In Sec. [5j infinity many different 2-adic wavelet orthonormal bases in £ 2 ((Q>2) 
are constructed. Namely, using Theorem 15.11 we construct wavelet functions 
ip( s \x), s £ N whose dilatations and shifts form 2-adic orthonormal wavelet 
bases in C 2 (Q 2 ). 

Since many p-adic models use pseudo-differential operators, in particular, 
fractional operator, these results on p-adic wavelets can be intensively used in 
applications. Moreover, p-adic wavelets can be used to construct solutions of 
linear and semi-linear pseudo-differential equations [5], [23] . 



2. p-ADlC DISTRIBUTIONS 

We recall some facts from the theory of p-adic distributions (generalized 
functions). Here and in what follows, we shall systematically use the notations 
and results from and [HI Ch.II]. Let N, Z, C be the sets of positive 
integers, integers, complex numbers, respectively, and No = {0} U N. Denote 
by Qp = Q P \ {0} the multiplicative group of the field Q p . 

The canonical form of a p-adic number x ^ is 

(2.1) x = p 7 (x + Xip + x 2 p 2 H ), 

where 7 = 7(2) £ Z, Xj — 0, 1, . . . , p — 1, x ^ 0, j — 0, 1, The series 

is convergent in the p-adic norm (ll.ip . and one has \x\ p = p -7 . By means of 
representation (12. ip . the fractional part {x} p of a number x £ Q p is defined as 
follows 

( 9 o\ f T \ = / °' if ^ or x = °' 

1 ' 1 ip \ p^x + x 1 p + x 2 p 2 + --- + x bhl p^- 1 ), if 7(x)<0. 

The function 

(2.3) Xp (&) = e 2 ^ x} " 

for every fixed ^ £ Q p is an additive character of the field Q p . 

According to [331 HI. 2.], any multiplicative character tt of the field Q p can 
be represented as 

Tt(x) d = Tc a (x) = \x\ p ~ l 'K l {x), x £ Q*, 



where ir(p) = p l ~ a and n\{x) is a normed multiplicative character such that 
7Ti(x) = 7Ti(|a;|pa:), 7Ti(p) = 7Ti(l) = 1, | ttx (a;) | = 1. We denote 7r = l^l" 1 . 

The space Q™ = Q p x • ■ ■ x Q p consists of points x = (xi, . . . , x n ), where 
Xj £ Q p , j — 1, 2 . . . , n, n > 2. The p-adic norm on QJJ is 

(2.4) \x\ p = max |x,-L, x £ Q™, 

l<j<n ^ 

where \xj\ p id defined by (II. ip . 
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Denote by B™(a) = {a; G Q™ : \x — a\ p < p 1 } the ball of radius p 1 with the 
center at a point a = (a 1; . . . , a n ) G Q™ and by 5™(a) = {x G Q™ : \x — a\ p = 
p 1 } = B y (a) \ 5"_ 1 (a) its boundary (sphere), 7 G Z. For a = we set 
£>" (0) = 5™ and S y (0) = S 1 ™. For the case n = 1 we will omit the upper index 
n. It is clear that 

(2.5) B^(a) = S 7 (a x ) x ■ • • x B 7 (a n ), 

where B 1 {aj) = {xj : |xj — | p < p 7 } C Q p is a disc of radius p 7 with the 
center at a point a,- G Q p , j = 1, 2 . . . , n. 

Any two balls in Q p either are disjoint or one contains the other. Every 
point of the ball is its center. 

According to [SJ 1. 3, Examples 1,2.], the disc B 7 is represented by the sum 
of p 7 " 7 disjoint discs By (a), 7' < 7: 

(2.6) B 7 = ByUU n B f (a), 

where a = and a = a_ r p~ r + a_ r+ ip~ T+1 + ■ — h a_y_ip _7 ' _1 are the centers 
of the discs By (a), r = 7, 7 — 1, 7 — 2, . . . , 7' + 1, < < p — 1, a_ r 7^ 0. 
In particular, the disc B is represented by the sum of p disjoint discs 

(2.7) J B = J B_ 1 UU^ 1 1 J B_ 1 (r), 

where £>_i(r) = {x & S : x = r} = r + pZ p , r — 1, . . . ,p — 1; £>_i = < 
P 1 } — P^pi an d 'S'o = {\ x \p — 1} = U^ 1 5_i(r). Here all the discs are disjoint. 
We call coverings (I2.6P and (I2.7P the canonical covering of the discs B Q and 
By, respectively. 

On Q p there exists the Haar measure, i.e., a positive measure dx invariant 
under shifts, d(x + a) = dx, and normalized by the equality <x dx = 1. 
The invariant measure dx on the field Q p is extended to an invariant measure 
d n x = dx\ ■ ■ ■ dx n on Q™ in the standard way. 

If / is an integrable function on Q p , then [TH Ch.II,§2.2], [33 IV]: 



dx = p 1 , 

(2.8) / f(x)dx = T. N 1= -oo! s J^)dx, 

I f(x) dx = J B / f(x)dx - S Bn _ 1 f{x)dx. 

A complex-valued function / defined on Q™ is called locally- constant if for 
any x G Q p there exists an integer l(x) G Z such that 

f(x + y) = f(x), yeBfc. 

Let £(Q p ) and V(Q p ) be the linear spaces of locally- const ant C-valued func- 
tions on Q" and locally-constant C-valued functions with compact supports 
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(so-called test functions), respectively [33, VI. 1., 2.]. If <p G £>(Q™), according 
to Lemma 1 from [33| VI. 1.], there exists I G Z, such that 

<p(x + y) = <p(x), yeB?, x G . 

The largest of such numbers I = l(ip) is called the parameter of constancy of 
the function <p. Let us denote by V l N (Qp) the finite-dimensional space of test 
functions from X>(Q") having supports in the ball and with parameters 
of constancy > / [33J VI. 2.]. The following embedding holds: Pjy(Q") C 
^(Qp)» N <N',l> V. Thus D(Q£) = limind 7V -,oolinimd^_ 00 P^(Q«). 
The space T>(Qp) is a complete locally convex vector space. 

According to [33J VI, (5. 2')], any function ip G V l N (Qp) is represented in the 
following form 

p n(N-l) 

(2.9) <p{x)= <p{<nMz-<n, ^q;, 

where Ai(x — c u ) are the characteristic functions of the disjoint balls Bi(c u ), 
and the points c u = (c£, . . . c£) G -Bjy do not depend on <p. 

Denote by P'(Qp) the set of a11 linear functional on £>(Q") [33J VI.3.]. 

Let us introduce in U(Q™) a canonical S-sequence S^x) = p nk Q(p k \x\ p ), and 
a canonical 1-sequence A^(x) = fi(p~ fc |x| p ), fceZ, a; G Q™, where 

(2.10) n(i) = {J; 05 If ]; 

Here A^(x) is the characteristic function of the ball B%. It is clear [33[ VI.3., 
VII.l.] that 6 k -»• 5, Jfe -»• oo in £>'(Q£) and A fe -»• 1, jfe -»• oo in £(<Q£). 
The Fourier transform of ip G X>(Q") is defined by the formula 

*M£)= / x P (£ • xV(x) cTx, ^gQ^, 

where Xp(£ ' x ) = XpC&^i) ' ' 'Xp{£,nXn) = e 2m ^i= 1 ^ jX ^ p ] £ • a; is the scalar 
product of vectors. 

The Fourier transform is a linear isomorphism T>(Q p ) into X>(Q"). Moreover, 
according to [321 Lemma A.], [33, 111,(3.2)], [Ml VII.2.], 

(2.11) ^x)eV l N (Q n P ) iff 

We define the Fourier transform F[f] of a distribution / G P'(Qp) by the 
relation [32 VII.3.]: 

(2.12) WU} = (/^M), v^gd(q;). 

Let A be a matrix and b G Q". Then for a distribution / G ^'(QL 1 ) the 
following relation holds [34, VII, (3.3)]: 

(2.13) F[f(Ax + 6)](0 = | det A\^ Xp ( - A~ x b ■ £)F\J(x)) {A~ x £), 
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where det A ^ 0. According to [34, IV, (3.1)], 

(2.14) F[A k }(x) = 6 k (x), kez, xeq;. 

In particular, F[f2(|^| p )](x) = f2(|x| p ). 

The convolution / * g for distributions f,g G 2?'(Qp) is defined (see [34], 
VII.l.]) as 

(2.15) (/ *g,<p)= lim (/(a;) x g(y), A k (x)(p(x + y)) 

k— >oo 

if the limit exists for all <p G T>(Q™), where f(x) x g(y) is the direct product 
of distributions. If for distributions /, g G V'(Qp) the convolution / * g exists 
then [3U VII,(5.4)] 

(2.16) F[/*#] = F[/]F[<7]. 

Definition 2.1. Let i\ a be a multiplicative character of the field Q p . A dis- 
tribution / G T>'(Qp) is called homogeneous of degree 7Tq, if for all if G 2}(Q") 
and t G Q* we have the relation 

(/. Kt' ' ' ' ' t) ) = 7Fa(t)|t| ^ / ' • • • ' 

i.e., f(tx) = f(txx,...,tx n ) = ir a (t)f(x), x = (xi,...,x n ) G Q™. A homoge- 
neous distribution of degree ir a (t) = |t| p _1 (a ^ 0) is called homogeneous of 
degree a — 1. 

3. The p-adic Lizorkin spaces 

Let us introduce the p-adic Lizorkin space of test functions 
$(Q^) = {0:0 = FM,^GVI/(Q«)}, 

where 

*(<Q£) = MO e ^ (Qp) = ^(0) = o}. 

Here #(Q£), $(Q") C D(Q£). The space $(Q") is called the p-adic Lizorkin 
space of test functions. The space $(Q p ) can be equipped with the topology 
of the space £>(Q") which makes $ a complete space. 
In view of (12. lip , the following lemma holds. 

Lemma 3.1. ( [3], [I]) (a) G $(Q£) iff (p E £>(Q") and 
(3.1) / 0(x)cTx = O. 

(b) G Pjv(Q^) n $(Q£), i.e., / B „ 0(0;)^ = 0, z# V = e 

Dlf (Q£) n tf(Q£), i.e., = o, e e S-jv- 

Unlike the classical Lizorkin space, any function ^(0 £ ^(Qp) is equal to 
zero not only at £ = but in a ball B n 3 0, as well. 

Let $'(Qp) denote the topological dual of the space $(Q p ). We call it the 
p-adic Lizorkin space of distributions. 
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By ty 1 - and $- L we denote the subspaces of functional in P'(Qp) orthogonal 
to #(<Q£) and $(Q"), respectively. Thus * x = {/ G W(Q%) '■ f = CS, C G C} 
and $-L = {/ G D'(Q£) : f — C, C E C}. 

Proposition 3.1. ( [3]) 

The space can be obtained from £>'(Qp) by "sifting out" constants. 

Thus two distributions in £>'(Qp) differing by a constant are indistinguishable 
as elements of $'(Q"). 

Similarly to (I2.12p . we define the Fourier transform of distributions / G 
$' X (Q") and g G tf' x (Q£) by the relations: 

OTLVO = </,*», v^g*(q-), 

(i^UH^M), V0G$(Qp. 

By definition, F[$(Q™)] = (Q£) and F[#(Q£)] = $(Q£), i.e., (Q give well 
defined objects. 

4. Construction of multiresolution analysis 

4.1. p-Adic multiresolution analysis. Denote the factor group Q p /Z p by 
I p , i.e. 

I p = {a = jT 7 (a + aip H h a 7 _ip 7_1 ) : 

(4.1) 7 G N; % = 0, 1, . . . ,p - 1; j = 0, 1, . . . , 7 - 1}. 

It is well known that Q p = B U U^ =1 S^, where 5 7 = {x G Q p : |x| p = p 1 }. 
In view of (12.11) . x G 5* 7 , 7 > 1 if and only if x = a;_ 7 p~ 7 + x_ 7+ ip~ 7+1 H — • + 
x-ip^ 1 + £, where £ G £>o- Since a;_ 7 p~ 7 + x_ 7+ ip~ 7+1 + ■ • • + X-i^r 1 G J p , we 
have a "natural" decomposition of Q p to a union of mutually disjoint discs: 

Q p = U aeIp B (a). 

So, /„ is a "natural" group of shifts for Q p . 

Definition 4.1. A collection of closed spaces Vj C £ 2 (Q P ), j G Z is called a 
multiresolution analysis (MRA) in £ 2 (Q P ) if the following axioms hold 

(a) Vj C V j+1 for all j G Z; 

(b) U ieZ ^- is dense in £ 2 (Q P ); 

(c) n jeZ Vj = {0}; 

(d) /(•) G Vj ^ /(p- 1 -) G V, +1 for all j G Z; 

(e) there a function G Vo such that the system <p(x — a), a G I p , form an 
orthonormal basis for Vq. 

The function from axiom (e) is called scaling or refinable. It follows im- 
mediately from axioms (d) and (e) that the functions p^ 2 (p(p~ : ' ■ —a), a G I p , 
form an orthonormal basis for Vj. 
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According to the standard scheme (see, e.g., [31;, §1.3]) for construction of 
MRA-based wavelets, for each j, we define a space Wj (wavelet space) as the 
orthogonal complement of Vj in Vj + i, i.e., 

(4.2) V j+1 = V 3 © Wj, j G Z, 
where Wj 1 j G Z. It is not difficult to see that 

(4.3) / G Wj -<=>• f{p 1 -) G Wj'+ij for all j'GZ 

and Wj _L Wk, j ^ k. Taking into account axioms (b) and (c), we obtain 

(4.4) ®j&Wj = £ 2 (Q P ) (orthogonal direct sum). 

If now we find a function ip G Wq such that the system ip(x — a), a G I p , form 
an orthonormal basis for Wq, then the system p^ 2, ip(p~ : ' ■ —a), a G I p , is an 
orthonormal basis for £ 2 (Q P ). Such a function ip is called a wavelet function 
and the basis is a wavelet basis. 

4.2. p-Adic refinement equation. Let be a refmable function for a MRA. 
As was mentioned above, the system p l l 2 cf){p~ l ■ —a), a G I p , is a basis for V±. 
It follows from axoim (a) that 

(4.5) = 2J cta^ip' 1 ■ -a), a a eC. 

a£l p 

We see that the function is a solution of a special kind of functional equation. 
Such equations are called refinement equations. Investigation of refinement 
equations and their solutions is the most difficult part of wavelet theory in real 
analysis. 

A natural way for construction of a MRA (see, e.g., (311 §1.2]) is the fol- 
lowing. We start with an appropriate function <p whose integer shifts form 
an orthonormal system, and set Vo = span {0 (x — a) : a G I p j and Vj = 

span{0(p~- 5 'a; — a) : a G J p }, j G Z. It is clear that axioms (d) and (e) of 
Definition 14.11 are fulfilled. 

Of course, not any such a function provides axiom (a). In the real setting, 
the relation Vq C V\ holds if and only if the refmable function satisfies a 
refinement equation. Situation is different in p-adics. Generally speaking, 
a refinement equation (14.51) does not imply the including property Vq C V±. 
Indeed, we need all the functions 0(- — b), b G I p , to belong to the space V±, 
i.e., the equalities 0(x — b) = J2 a ei P a a,b4>(p~ lx ~ a ) should be fulfilled for all 
b G I p . Since p~ 1 b + a is not in I p in general, we can not state that refinement 
equation (14.51) implies 0(x — b) = J2 a ^i a a,b4>(p~ lx ~ P~ x b — a) G V\ for all 
b G I p . 

The refinement equation reflects some "self-similarity" . The structure of the 
space Q p has a natural "self-similarity" property which is given by formulas 
(12. 6p . (12.71) . By (12.71) . the characteristic function A (x) = f^(|x| p ) of the unit 
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disc B is represented as a sum of p characteristic functions of the disjoint 
discs B_i(r), r = 0, 1, . . . , p — 1, i.e., 

(4.6) A 0) = VA (-x--), xgQ p . 

r=0 ^ ^ 

Thus, in p-adics, we have a natural refinement equation (14.51) : 

P_1 1 r 

(4.7) ^^Vrf-^--), xeQ p , 

whose solution is <f)(x) = A (x) = This equation is an analog of the 

refinement equation generating Haar MRA in real analysis. 

4.3. Construction of 2-adic Haar multiresolution analysis. Now, using 
the refinement equation (14. 7p for p = 2 

(4.8) = + - |), x£®2, 

and its solution, the refinable function <p(x) = Aq(x) = fH|:c| 2 ), we construct 
2-adic multiresolution analysis. 
Set 

(4.9) Vq = span{0(x — a) : a E J 2 }, 

(4.10) Vj = span{0(2-^x- a) : a G i 2 }, jeZ. 

It is clear that axioms (d) and (e) of Definition 14.11 are fulfilled and the system 
2^ 2 0(2-^' ■ -a), a G I p is an orthonormal basis for V^-, j G Z. 

Note that the characteristic function of the unit disc fM|a;| 2 j has a wonderful 
feature: Q(\ ■ +£| 2 ) = ^(1 ' b), for all £ G Z 2 because the p-adic norm is non- 
Archimedean. In particular, 0(| • ±1| 2 ) = fl(\ ■ | 2 ), i.e., 

(4.11) <j>(x±l)=<p(x), VxgQ 2 . 

Thus (f> is periodic with the period 1. 

In view of this fact, taking into account that 2 _1 6 + a ( mod 1) is in 7 2 , for 
all a, b G 7 2 , it follows from the refinement equation ( 14.81) that Vo C Vi. By 
( 14.101) . this yields axiom (a). 

Due to the refinement equation (14.81) . we obtain that Vj C Vjr+i, i.e., the 
axiom (a) from Definition 14.11 holds. 

Lemma 4.1. The axiom (b) of Definition 14.11 holds, i.e., Uj^zVj = £ 2 (Q 2 ). 

Proof. According to (12.91) . any function if G T>(Q 2 ) belongs to one of the spaces 
D l N (Q.2), and consequently, is represented in the form 

p N-l 

(4.12) <p(x) = ^#)A,(x-0, x G Q 2 , 
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where Aj(- — c v ) are the characteristic functions of the mutually disjoint discs 
Bi{c v ) C Q 2 , c v e B N , v = l,2,...p N ~ l ; I = l(<p), N = N(<p). Since 
Ai(x — c u ) = fl(p \x — c u \ p ) = Q(\p l x — p l c u \ p ) and any number p l c v can be 
represented in the form p l c v = a u + b v , where a? G I 2 , b u G Z 2 , we have 
Aj(x — c") = A^(x — a y ). Thus any function y9 G P(Q2) can De represented in 
the form 

p N-l 

(4.13) = ^2a„Ai(x-a v ), x G Q 2 , a u E I 2 , a u G C. 

i/=l 

Consequently, on the basis of (I4.10p . yj(a;) G Thus any test function ip 
belongs to one of the space Vj, where j = j(<f). 

Since the space T>(Q 2 ) is dense in £ 2 (Q 2 ) [34, VI. 2], approximating any 
function from £ 2 (Q 2 ) by test functions (14.131) . we prove our assertion. □ 

Lemma 4.2. The axiom (c) of Definition 14.11 holds, i.e., Hj^zVj = {0}. 

Proof. Suppose that Dj^zVj 7^ {0}- Then there exists a function / G Vj for all 
j G Z. Hence, due to (HlDjl . /(re) = Eae/ 2 c ia0(2 _i z - a) for all j G Z. 

Let x = 2- N (x +x 1 2+x 2 2 2 + - ■■). Since 2^x = 2~ N ^{x +x 1 2+x 2 2 2 + - ■ ■ ), 
for all j < —N, we have 2~^x G Z 2 , and, consequently, |2~ J x — a\ 2 > 1 for 
all a G J 2 , a 7^ 0. Thus (j){2~^x — a) = for all j < —V and a G J 2 , a 7^ 0. 
Since [ 2 — J a? 1 2 < 1, we have f(x) = Cjo for all j < —N. Similarly, for another 
x' = 2- N '(x' + x\2 + x' 2 2 2 + •••), we have f(x') = c ro for all j < —N'. This 
yields that f(x) = f(x'). Consequently, f(x) = C, where C is a constant. 
However, if C 7^ 0, / ^ £ 2 (Q 2 ). Thus, C = and the proof of the theorem is 
complete. □ 

According to the above scheme, we introduce the space Wo as the orthogonal 
complement of Vo in V\. 
Set 



(4.14) 




Lemma 4.3. The shift system ip(°\x — a) , a G I 2 , is an orthonormal basis of 
the space W . 

Proof. Let us prove that W _L V Q . It follows from (Oji . (Oil) that 



(^ (o) (x-a),0(x-6)) = / ^ (0) (a; - a)0(x - 6) da; 




for all a, b G I 2 . Let a ^ b. Since it is impossible & ^a + l, taking 

into account that the functions 2 1//2 0(2~ 1 • — c), c £ I 2 are orthonormal, we 
obtain — a), <p(x — b)) =0. If a = b, again due to the orthonormality 
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of the system 2 1/,2 0(2 1 • — c), c G I2, taking into account that |, | + | G I2, 
we have 




Thus, ip(°\x + a) _L 0(x + 6) for all a, b G I2. 

The refinement equation ( 14. 8f) and relation (14.141) imply that 



Since {2 1 / 2 0(2 _1 x — a) : a £ I 2 } is a basis for V\, we have Vi = Vo © Wo, i-e., 
(J42D holds. □ 

Thus we prove that the collection {Vj : j G Z} is a MRA in £ 2 (Q2) and 
the function ^ (0) defined by gUj) is a wavelet function. This MRA is a 2- 
adic analog of the real Haar MRA and the wavelet basis generated by ip^ is 
an analog of real Haar wavelet basis. But in contrast to the real setting, the 
refinable function <fi generating our Haar MRA is periodic with the period 1 
(see (14.111) ). which never holds for real refinable functions. It will be shown 
bellow that due of this specific property of 0, there exist infinity many different 
orthonormal wavelet bases in the same Haar MRA (see Sec. [5]). 

Due to (JZ2]), flUD, the function ^ (0) can be rewritten in the form ^°'(x) = 
X2^ 1 x)VL{\x\2) and the Haar wavelet basis is 

(4.15) =2-~ l/2 X2 {2- 1 (Vx-a))tt(\T 1 x-a\ 2 ), x G Q 2 , 7 e Z, a G J 2 . 
It is clear that 

(4.16) / ^(z)cb = 0, 

and, according to Lemma [3.11 ip$(x) belongs to the Lyzorkin space $(<Q>2)- 

Remark 4.1. The Haar wavelet basis (14. 15|) coincides with Kozyrev's wavelet 
basis (11.21) for the case p = 2. In present paper we restrict ourself by con- 
structing the Haar wavelets only for p = 2. Since Haar refinement equation 
(14. 7p was presented for all p, a similar construction may be easily realized in 
the general case. Moreover, it is not difficult to see that Kozytev's wavelet 
function 9j(x) from ( jl.2p can be expressed in terms of the refinable function 
4>(x) as 

(4.17) 9j{x) = Xp (p- 1 Jx)n(\x\ p ) =p- 1/2 J2 h r<p(-x--), ^Q P , 
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where h r = p l / 2 e 2m ^- p r — 0, 1, . . . ,p — 1, j = 1, 2, . . . , p — 1. 

Remark 4.2. In view of periodicity (14.111) of the refinable function 0, one can 
use shifts ip(°'(- + a), a G h, instead of shifts — a), a G -?2- 

Now we show that there is another function ij)( l >(x) whose shifts form an 
orthonormal basis in Wq. Indeed, taking into account (14.111) . we have 




4) {1) {x-a) 

1 / fx a\ fx a 1\ fx a 1\ fx a 1 \\ 

1 / , fx a\ ,fx a 1\ , fx a 1 1\ , /x a 1 \\ 
(420) =7ir(2-2)-H2-2-2)-<2-2-25-2) + H2-2-25)j' 

Since the system of functions {<f)(2^ 1 x — a) : a G J2} is orthonormal, in 
view of (14.111) . formulas (14 . 1 8 j) — (14 . 20 j) imply that the function ip^(x) and the 
function ip^(x — a) are orthonormal, whenever a G I2, a 7^ 0, ~. Here we 
take into account that all shifts (up to mod 1) of refinable function in (14.181) . 
(H20D are distinct. 

Similarly, by (I4TTHD . (Q5|) . we have 

(^ {1) (x),^ 1) (x + 2- 1 )) = [ ^ 1 \x)^ 1 \x + 2~ l )dx 



and 




14 



V. M. SHELKOVICH AND M. SKOPINA 



Thus all shifts of ip^ 1 ' are orthonormal. 

It is clear that the functions (14.181) and (I4.19P can be rewritten in the form 



(4.21) 




It follows that 




Since the system vp(°\- — a), a G I2, forms an orthonormal basis for Wo, the 
system ^>W(- — a), a G I2, is another orthonormal basis for Wo. 

So, we showed that a wavelet basis generated by the Haar MRA is not 
unique. 

5. Description of 2-adic Haar bases 

5.1. Complex wavelets. Using the fact that all dilatations and shifts (x — > 
2 7 x + a, a G I2) of the Haar wavelet function form a orthonormal basis in 
£ 2 (Q2), we show that there exist infinitely many wavelet functions ^ > s G N 
in W . 

In what follows, we shall write the 2-adic number a = 2~ s [ao + a\2 + • • • + 
a s _i2 s_1 ) G I2, dj = 0, 1, j = 0, 1, . . . , s — 1 briefly as a rational number 
a — ||, where m = a + d\1 + • ■ ■ + a s _i2 s_1 . 

Since the characteristic function of the unit disc <p(x) = A (a;) = 0(|x|2) is 
periodic with the period £ G So, the wavelet function ifj°(x) has the following 
evident and important property: 

(5.1) ?P {0) (x + = -iP {0) (x), £eSo. 

Here £ = 1 + £2 + £ 2 2 2 + • • • , where ^ = 0, 1; j G N. 

Before we prove a general result, we consider the simplest particular case. 
Consider the function 

(5.2) ip {1 \x) = ao^J {0) (x) +a 1 ip {0) (x + %«i6C, 

and solve the problem when all shifts of this function generates an orthonormal 
basis ip^(x + a), a G h in Wo- 

Taking into account orthonormality of the system ip^(x + a), a G I2 and 
relation (15. ip . we can see that the function ip^(x) and the functions ip^(x + a) 
are orthonormal for all a G I2, a ^ 0, |. Thus, in view of (15.11) . the system 
of functions ip^ l '{x + a), a G -Z2 is orthonormal if and only if the system of 
functions (15. 2p and 

(5.3) (x + -) = -ai^ (0) (x) + « ^ (0) (z + -) 
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is orthonormal. Hence, we have |ao| 2 + |cti| 2 = 1- In other words, the matrix 

is unitary. Thus, the function (15.21) . where |«o| 2 + |«i| 2 = 1 is the wavelet 
function. It is clear that the wavelet function (14.2ip is a particular case of the 
wavelet function (15.21) . 

Consequently, all dilatations and shifts of ip^(x) form 2-adic orthonormal 
wavelet basis in £ 2 (Q 2 ). 

Now we will prove a general theorem. 



Theorem 5.1. Let s = 1, 2, . . . . The function 

2 S -1 

(5.4) 

fc=0 

is the wavelet function (whose dilatations and shifts form 2-adic orthonormal 
wavelet basis in £ 2 (Q 2 ) ) if and only if 



^' ) (s) = X>^< D) (* + £ 



2 S -1 



(5.5) a k = 2" 

lk e C, |7 fc | = 1. 



, fc = 0,l,2,...,2' 



r=0 



Proof. Suppose that ip^ s '(x), s > 1 is given by formula (15 Ah . Since the system 
+ a), a e I2 is orthonormal (see Subsec. 14. 3[) and in view of relation 
(15. ip . it is easy to see that ip^ and + a) are orthonormal for any a E I2, 

a 7^ Jh k = 0, 1, . . . 2 s — 1. Thus the system of functions ip^ s \x + a), a e I 2 is 
orthonormal if and only if the system of functions, consisting of the function 
(15 .4p and its shifts, i.e., 

^ {s) {x + ^) = -a 2 s-r^\x)-a2s- r+ ^\x+^ « 2s _^(°)(x+^) 



(5.6) 



+ a ^ (0) ( x + — ) H h a^-r-iV'^ ( & + 



(0) 



r = 0, 1, . . . , 2 s — 1 is orthonormal. 

Set S(°> = {^ (0) (- + |) : k = 0,1,...,2 S -1} T , H« = {V> (s) + = * 
0, 1, . . . , 2 s - 1} T . In view of (I5~4l . = DE^°\ where 

/ 



(5.7) 



D 



V 



Q'o 


a 1 


«2 


Oi2 s ~2 


a 2 s_i 


— OL 2 *-l 


Of, 






Q>2 s -2 


—a<2< i -2 


— 0:2^-1 


a . . 


ct2 s -4 




—OL2 




— Gi4 . . 


a 


Oil 




— 0:2 


-a 3 . . 


. — a^-i 


OiQ 



Thus the system 5^ is orthonormal if and only if the matrix D is unitary. 
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Let u = (a , oti, . . . , a 2 s -i) T be a vector and 





f° 


. 


. 





-1 


\ 




1 


. 


. 










A = 





1 . 


. 



















. 


. 1 












Vo 


. 


. 


1 





/ 



be a 2 s x 2 s matrix. It is easy to see that 

A r u = (-a 2 s- r , ~a 2 s_ r+1 , • • • , —a 2 s-i, «o, a 2 s- r -i) 1 ', 

r = 1, 2, . . . , 2 s — 1. Thus D = (u, Au, . . . , A^^u) 1 . It is significant that 
A 2S u = —u. Consequently, in order to describe all matrixes D (or in other 
words, all vectors u), we should find all vectors u = (a , a±, . . . , a 2 ^i) T such 
that the system {A r u : r = 0, 1, 2, . . . , 2 s — 1} is orthonormal. 

In view of the fact that the system ip^(x + a), a G I 2 forms an orthonormal 
basis in W , it is easy to see that the vector u = (1,0,..., 0, 0) T is one of 
mentioned above vectors u. That is the system composed of vectors uo and 
A r u = (So r , Sir, ... , ^2 a -2r, ^2 a -irY ', r = 1; 2, ... , 2 s — 1, is orthonormal, where 
5i r is the Kronecker symbol. 

Let us prove that the vector u = («o, «i, • • • , «2 s -i) T already mentioned 
above such that A r u, r = 0, 1, 2, . . . , 2 s — 1 is orthonormal, can be expressed 
by the formula u = Buq if and only if B is a unitary matrix such that AB = 
BA. Indeed, let u = Bu , where B is a unitary matrix such that AB = 
BA. Then A r u = BA r u , r = 0, 1, 2, . . . , 2 s - 1. Since the system A r u , 
r = 0, 1, 2, . . . , 2 s — 1 is orthonormal and the matrix B is unitary, the vectors 
A r u, r = 0, 1, 2, . . . , 2 s — 1 are orthonormal. Conversely, if the system A r u, 
r = 0, 1, 2, . . . , 2 s — 1 is orthonormal, taking into account that the system A r u , 
r = 0, 1, 2, . . . , 2 s — 1 is orthonormal, we conclude that there exists a unitary 
matrix B such that A r u = B(A t uq), r = 0, 1, 2, . . . , 2 s — 1. Since A 2 ' 'u = —u, 
A 2S u = —u , we have an additional relation A 2 " u = BA 2S u$. It follows from 
the above relations that {AB - BA)(A r u ) = 0, r = 0, 1, 2, . . . , 2 s - 1. Since 
the vectors A r u , r = 0, 1, 2, . . . , 2 s — 1 form a basis in the 2 s -dimensional 
space, we conclude that AB = BA. 

Thus we have D = (Buq, BAuo, . . . , BA 2S ~ 1 uq) . 

It is clear that the eigenvalues of A and the corresponding normalized eigen- 
vectors are 

(5.8) A r = -e* 7 ^, 

and v r = ((f r )i, . . . , (v r ) 2 s^ T , respectively, where 



(5.9) 



{v r ) t = 2-*/ 2 (-l)' e — I = 0, 1, 2, . . . , 2 s - 1, 
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r = 0, 1, 2, . . . , 2 s — 1. As is well known, the matrix A can be represented as 
A = CAC' 1 , where 

/ A . . . \ 
Ai ... 



.4 



/ 



\ ... A 2s _i 

is a diagonal matrix, C = (v ,Vi, . . . ,v 2 ^-i). Since C is a unitary matrix, the 
matrix B = CBC~ l is unitary if and only if B is unitary. On the other hand, 
AB = BA if and only if AB = BA. Moreover, since according to (15.81) A& ^ A;, 
whenever k ^ I, all unitary matrix B such that AB = BA, are given by 

\ 

7! ... 

B 



V ... 72s _! / 



where 7^ G 



iTfcl 



1. Hence, all unitary matrix B such that AB = BA, are 
given by B = CBC~ X , where B is the above diagonal matrix. 
By using formula ( 15.91) . one can calculate 



a k = (Bu ) k = {CBC l u Q ) k = ^2 lr{v r . 



k{Vr)0 



r=0 



_2r + l. 



k = 0,1,2, 



2 s - 1, 



r=0 



where 7^ G C, |7fe| = 1. Thus (I5.5P holds. 

Taking into account that S*- -* = D _1 H( S \ we conclude that if we define 
tjj^(x) by formula (15 .4p . where is given by (15. 5p . k = 0,1,2, ... ,2 s — 1, 



then the system of functions 



')(. - 



a £ a £ 1 2} is orthonormal and 



forms the orthonormal basis in Wq. 

Consequently, all dilatations and shifts of the function (15. 4p form 2-adic 
orthonormal wavelet basis in £ 2 (Q 2 ). □ 



It is clear that L ip$(x) dx = 0, and in view of Lemma IXTj tjj$(x) belongs 



(«), 



to the Lizorkin space G $( 



5.2. Real wavelets. Using formulas (15.51) . one can extract all real wavelet 
functions (15. 4p . 

Let s = 1. According to ([ED, Q , 



(5.10) 



(2) = cos #V> (0) (V) + sin#^ (0) ( x H — 



is the rea/ wavelet function. 
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Let s = 2. Set 7 r = e i0r , r = 0, 1, 2, . . . , 2 s - 1. Then (Ell) imply that the 
wavelet function ip^(x) is real if and only if 

sin 6*i + sin 2 + sin 3 + sin 4 = 0, 

cos 6*i — cos 6 2 + cos #3 — cos #4 = 0, 

sin 6\ — sin 2 — sin 3 + sin 4 = 
cos 6*i + cos 6*2 — cos #3 — cos 6*4, 

sin 9% — sin 2 — sin 83 + sin 4 = 
— (cos#i + cos 6 2 — cos 6*3 — cos 6*4). 

The last relations are equivalent to the system 

sin 6*i = —sin #4, cos 61 = cos 6*4, 
sin 6*2 = — sin0 3 , cos 6 2 = cos0 3 . 

Thus for s = 2 the real wavelet functions (15.41) is represented as 
^ l \x) = i(cos#i + cos0 2 )V (o) (x) 

+ ^=(cos#i - cos 9 2 + sin 0i + sm0 2 )^ (o) (x + 
+^ (sin 0i -sin0 2 )^ (o) (^+2) 

(5.11) + — (cos 0i - cos 2 - sin 0i - sin0 2 )^ (o) (x + — + 

In particular, for the special cases 0i = 2 = 0, 0i = —02 = 0, 0i = 02+f = 0, 
we obtain one-parameter families of the real wavelet functions 

+ sin 07/>(°) ( x + \ ) , 
(x) = cos 0^ (o) (x) + ^= sin 0</> (o) (x + ^) 

(5.12) _^ sin ^(o)^ + ^ + I^ 

^«( x ) = i(cos0-sin0)^(°)(x) + i ^(cos0 + sin0)^°)(x + ^) 

-i(cos0-sin0)^ {o) (x + i), 

respectively. 
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